Mechanism of thermally activated c-axis dissipation in layered High-T c 

superconductors at high fields 
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We propose a simple model which explains experimental behavior of c-axis resistivity in layered 
High-T c superconductors at high fields in a limited temperature range, ft is generally accepted that 
the in-plane dissipation at low temperatures is caused by small concentration of mobile pancake 
vortices whose diffusive motion is thermally activated. We demonstrate that in such situation a 
finite conductivity appears also in c-direction due to the phase slips between the planes caused by 
the mobile pancakes. The model gives universal relation between the components of conductivity 
which is in good agreement with experimental data. 
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A well known property of layered High-T c supercon- 
ductors (e.g., Bi2SrCaiCuiO x (BiSCCO)) is an exis- 
tence of extended low-temperature tails in the tempera- 
ture dependencies of the both in-plane and out-of- 
plane |^-|^] components of resistivity for magnetic field 
applied along the c-axis. An important experimental fact 
is that the thermally activated behavior for both p a \, and 
p c is characterized by the same value of activation energy 
(see Ref. and Fig. [|) . This indicates that dissipations 
in both directions have essentially the same origin. 

The origin of the tail for the in-plane resistivity com- 
ponent is qualitatively understood. It is theoretically 
expected that at fields larger than the crossover field 
B cr w $o/(7s) 2 and temperatures larger than the de- 
coupling temperature Td ec the system can be treated as 
a set of almost independent two-dimensional (2D) vor- 
tex lattices in the layers || (here 7 is the anisotropy 
of the London penetration depth and s is the interlayer 
spacing). The in-plane dissipation can therefore be at- 
tributed to the thermally activated motion of point de- 
fects in the 2D lattices (weak disorder) or weakly pinned 
mobile pancake vortices (strong disorder) 0. For a 2D 
superconductor with weak pinning a crossover between 
free flux flow and thermally activated motion should take 
place at the melting temperature of the vortex lattice 
T m 1=3 0.014s$q/(47tA) 2 as observed in the thin films of 
a — Nb^Ge Q and in numerical simulations |lC| ]. This 
seems not to be the case for available samples of the 
BiSCCO compound where crossover to the thermally ac- 
tivated behavior at high fields (> 1 T) takes place at 
temperatures 50-70 K which is significantly larger than 
the estimate for the melting temperature T m w 13 AT ob- 
tained using experimental values of the London pene- 
tration depth A a b s» 1800A and s = 15A. This indicates 
that in available BiSCCO single crystals pinning is strong 
enough to destroy the melting transition in the field range 
1-10 T and the crossover temperature between free flux 
flow and thermally activated regime ( 50-70 K) has to 
be interpreted as the depinning temperature of a single 



pancake vortex renormalized by the intervortex interac- 
tions 0. In this case a very natural assumption is that 
the lattice is strongly pinned and the in-plane dissipation 
is caused by the motion of mobile pancake vortices with 
the small concentration n and diffusion constant D. In- 
plane conductivity within this model can be estimated as 
follows. In-plane transport current j drives mobile pan- 
cakes with velocity v = [is(<f>o/c)j. Here p, is the average 
pancake mobility, which is connected with the diffusion 
constant D by the Einstein relation, D = pT. Pancake 
motion produces electric field E — ($o/c)un. Therefore 
the conductivity in afe-planc o- a (, in this model is given by 



Tab 



c 



T 
snD 



(1) 



Both n and D presumably have Arrhenius-type temper- 
ature dependencies but their temperature dependence as 
well as detailed mechanism of in-plane motion have no 
influence for our further consideration. 

In the region where the conductivity in the c-direction 
is much larger than the normal conductivity it has the 
Josephson origin and is caused by phase slips between 
neighboring layers. It was assumed j|J^] that these phase 
slips are similar to the phase slips in a single Josephson 
junction with the typical area &o/B and therefore the 
Ambegaokar-Halperin formula for the resistivity of small 
Josephson junction can be used to describe the temper- 
ature and field dependencies of p c . Unfortunately this 
physically transparent assumption does not have micro- 
scopic support. 

Diffusive motion of pancake vortices at high enough 
temperatures leads to the time variations of the phase dif- 
ference between the layers. Such pancake-induced phase 
slips provides a candidate mechanism for the c-axis dissi- 
pation which was not considered so far. In this Letter we 
obtain c-axis conductivity due to such mechanism. The 
key assumption of our model is that the phase slips are 
caused the diffusive motion of a small amount of pancake 
vortices independently in different layers. The model is 
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applicable above the decoupling temperature where the 
global superconducting coherence is broken. This means 
that the model can be used to describe dissipation in a 
limited temperature range (~ 40-70 K). 

The starting point of our calculation is the Kubo for- 
mula which relates a c with the correlation function of the 
Josephson current: 
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Here jj is the Josephson current, 8<p — <j>2 ~ 4>i — -§f A z is 
the Gauge invariant phase difference between neighbor- 
ing layers; the phases <^i,2(r, t) are mainly determined 
by the vortex positions Ri^i in the layers, 0i,2(r, t) — 
J2i ( l ) v{^ — 'R-i,2i{t)) where 4> v {r) = atan(y/x) is the phase 
distribution around the core of a single vortex. We calcu- 
late a c in the lowest order with respect to the interplane 
Josephson coupling. Within this approximation one can 
neglect correlations between Josephson currents in differ- 
ent layers, which are already excluded in Eq. (|J). We can 
also neglect interplane phase correlations in the sine-sine 
correlation function in Eq. (|2|) . This means that averages 
like (cos 54>) assumed to be zero and 

(sin <ty(0,0) sin 5<f>(r,t)) « (l/2)Re[5(r,t)], 
with S(r, t) = (exp i (5cj>(r, t) - (50(0, 0))) . 

As follows from Eq. (|J) the value of the out-of-plane con- 
ductivity is determined by how fast the phase correlations 
decrease in space and time. The origins of decay of the 
phase correlations in space and time are very different. 
The decay in space is determined by the randomness in 
the vortex arrangement while the decay in time is deter- 
mine by the vortex motion. In our model the static phase 
difference (50(0, 0) — 5<f)(r, 0) is established by the almost 
all rigidly pinned vortices |IlJ while dynamic behavior 
is caused by the small amount of mobile vortices. In 
such a situation the correlation function can be splitted 
into the separately averaging static and dynamic parts, 
S(r,t) w S(r)S(t) and Eq. (§) becomes, 
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with 



S(r) = (exp [i ((50(0, 0) 
S(t) = (exp[i ((50(0,0) 



<50(r,O))]), 
6<t>(0,t))]) . 



In absence of correlations between vortex positions in 
neighbor layers the static part is expected to decay at 
distances of the order of the average intervortex spacing 
do- This assumption can be supported by direct calcu- 
lation for the case of randomly placed pancakes, which 
gives O 



S(r) » exp (-Tm v r 2 h\(r max /r)) 



(4) 



with n v = B/<£>q. We assume that the phase at given 
point is not sensitive to pancake displacements at dis- 
tances larger than some phase coherence distance r max , 
which provides upper cutoff to the logarithmically di- 
verging integrals. The static phase correlation function 
for randomly placed pancakes Eq. (0) actually has the 
fastest possible descend. Correlations in the pancake ar- 
rangements always increase the spatial range of S(r) . De- 
cay of the dynamic part is determined by the diffusion of 
mobile pancakes and is given by 

S(t) = exp {-2imDt ln(r 2 max /Dt)) (5) 

Substituting Eqs. (@) and (§) into Eq.|| we obtain |]| 
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Up to logarithmic factors and numerical constant this 
result immediately follows from Eq. (En assuming that 
the interlayer phase difference is correlated at distances 
of the order of the average vortex spacing and at times 
of the order of the typical "phase slip" time = 1/nD. 
Using the relation jj = and Eq. (Q) the last 

equation can be transformed to 
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with C — 
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ln(nrj^ ax ) ln(r max /o ) ' 



Therefore within our 



model the ratio of conductivities depends only upon ma- 
terial parameters, temperature, and field, and does not 
depend upon the details of the mechanism, responsible 
for dissipation. Eq. (Q) represents the main result of our 
Letter. 

The model can be verified independently by measur- 
ing the frequency dependence of <j c . The frequency de- 
pendence can be taken into account by adding the fac- 
tor cos(bjt) under the integral in Eq.(Q). This gives 
Lorentzian frequency dependence a c {oS): 
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The typical "phase slip" frequency 
Lo ps = 2nnDhi(nr 2 nax ) can be directly related to the in- 
plane component of resistivity p a b = 1 /<?ab using Eq. (p]) , 
which gives uj ps w 2 10 s [1/s] • T[K] • p a b[^ ■ cm]. The di- 
rect correlation between typical frequency in cr c (w) and 
p a b provides additional experimental check of the theory. 

Let us discuss the limitations of our model. At low 
enough temperature slowing down of the pancake motion 
in the layers leads to increase of the interlayer correla- 
tions so that the diffusion of pancakes in different layers 
can not be considered independent any more. Quanti- 
tatively the influence of the Josephson coupling on the 
diffusive motion is characterized by value of the Joseph- 
son force acting on pancake. This force is determined by 



2 



the variation of Josephson energy under small pancake 
displacement and can be calculated as 



/,7a (R, t) — — Ej 



^(r-Rft)) 

dr a 



sin 60(r, i) (9) 



Neglecting again the interplane phase correlation we ob- 
tain the correlation function of this force 
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irEj ln(r max /fl) 
n v ln(r max /a) 



S(t) (10) 



As one can see from this expression the Josephson force 
has typical amplitude fj ~ aEj and is correlated at 
times ~ t ps and at large distances of the order of r max . 
This means that mobile pancakes within the correlation 
area ~ r 2 nax will drift under this force with approxi- 
mately the same velocity v rj fifj during the "phase slip" 
time t ps . This drift motion will produce an extra phase 
change 6(f> w nr max fifjt 
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naxfj/T. One can neglect 
the Josephson coupling and treat pancake motion in dif- 
ferent layers independently when this extra phase shift 
is small, which gives the condition r max fj <C T. More 
precise estimate requires evaluation of the in-plane phase 
coherence distance r m ax- This is beyond the scope of the 
simple model proposed in this Letter. 

Several mechanisms limit applicability of this model at 
high temperatures. Quasistatic phase fluctuations with 
the amplitude x = ((^0) 2 ) /2 lead to additional sup- 
pression of the Josephson energy, which can be described 
by the Debye- Waller factor E e J^ « Ejexp(-x), and 
reduces the ratio <7 c /cab as compared to Eq.(|^). Two 
contributions to \ come from thermal oscillations of the 
pinned pancakes (xi) and from the spontaneous inter- 
plane phase fluctuations not related to pancake motion 
(X2), X — Xi + Xi- We estimate this contributions as 
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where oll ~ s<f>oj c /(c£) is the Labush parameter, j c is 
the critical current, and J = s<I>g/(7r(47rA) 2 ) is the in- 
plane phase stiffness. Taking j c = 10 6 A/cm 2 , T = 60K, 
B = IT, lnfr^jj^ai/T) = 5 we obtain the estimate xi ~ 
0.15. Using parameters of BiSCCO and taking the loga- 
rithmic factor in xi as 5 we obtain xi ~ 4[K]/(T C — T). 
Therefore the spontaneous interplane phase fluctuations 
can be ignored at temperatures ~ 10 K below T c . 

Another important factor, which can not be ignored 
at high temperatures, is the quasiparticle conductivity. 
This conductivity has to be added to the phase slip con- 
ductivity (^) and therefore it enlarges the conductivity 
ratio in comparison with prediction of Eq. (^) . 

To verify the prediction of the model we plot in Fig. |^ 
temperature dependencies of the conductivity ratio at 
different fields (the data || are the same as in Fig. ([!])). 
One can see that within the temperature range 50-70 K 



for fields 0.5-2 T a general trend is in agreement with 
Eq. (Q), i.e. the ratio decreases with temperature and 
field. Taking the conductivity anisotropy at 60 K and 1 T 
(« 2.7- 10 -6 ) and using the relation between the Joseph- 
son coupling energy and the Ginzburg-Landau anisotropy 
ratio 7, Ej = <i> 2 l /7r(47rA7) 2 s, we estimate from Eq. (0) 
the value of 7 as 370 which is within the range reported 
in the literature jl4j . An abrupt increase of the ratio at 
low temperatures indicates probably the emergence of the 
global coherence in the system and the glass transition. 

To verify quantatively the validity of Eq. (^) we need 
to know the temperature dependence of the Josephson 
coupling energy. This dependence is determined by the 
temperature dependencies of the tunneling integral t and 
superconducting order parameter 'J, Ej oc i|^| 2 . Usually 
the first dependence is neglected. However, as recently 
was proposed by Bulaevskii et al fl5|| , the semiconduct- 
ing temperature dependence of the normal-state c-axis 
conductivity a cn at low temperatures can be ascribed to 
the temperature dependence of t. We therefore assume 
that t and a cn have the same temperature dependence. 
The temperature dependence of a cn have been recently 
measured by Cho et al || for the temperatures down to 
sa 55 K using high magnetic fields (up to 18 T). This de- 
pendence in the interval 55-90 K can be very well fitted 
by the exponential formula <r cn cx exp(r/29[i4T]). There- 
fore to verify Eq. (0) we try to fit the conductivity ratios 
by the formula 

a c /a ab (B, T) = A(B) exp(2T/29)(l - T/T C ) 2 T~ 2 (11) 

with A(B) being the fitting parameter. The factor 
1 — T/T c accounts for temperature dependence of |^| 2 . 
As one can see this formula gives very good fit to the data 
in the temperature range 45-70 K. However the field de- 
pendence of the fitting parameter A(B) is weaker than 
predicted by Eq. (@), A{B) cx S" ' 52 . The most prob- 
able origin of this discrepancy is that increasing of the 
field reduces randomness in vortex arrangement within 
the layers and, as a consequence, the static phase corre- 
lations decay in space slower than is predicted by Eq. (Q). 

In conclusion, we propose simple model, which pro- 
vides the mechanism for thermally activated dissipation 
along the c-axis in the vortex state of disordered layered 
superconductors. The model gives the universal relation 
between the conductivity components. 
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FIGURE CAPTIONS 

FIG. 1. Temperature dependencies of the resistivity com- 
ponents p ab and p c for BiSCCO for fields 0.5, 1,2, and 5 T. 
Data are taken from Ref. Q. Note that in spite of huge dif- 
ference in absolute values, p a b and p c have almost identical 
temperature dependencies. 

FIG. 2. Temperature dependencies of the conductivity ra- 
tio for the data from Ref. Solid lines give the fits by 
the formula ([H]), which were made to verify the prediction 
of Eq.(Q). Table gives the values of parameter A at different 
fields obtained from the fits. 
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